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Abstract 


In t his t hesis, a comprehensive study on Dielectric itesemutor operat- 
ing in TM 0 is inode, kept in various types of environments, has been done. 
The analysis presented here, is applicable to Isolated DR, DR post, DR in a 
cavity, DR in cylindrical metal waveguide below cut-off, suspended substrate 
etc. Various parameters, required for the design of microwave circuits includ- 
ing DR for millimetcric wave applications, have been studied and analytical 
formulations have been developed for them. The resonance frequency of a 
cylindrical DR for TMqu mode, is calculated using Dielectric Waveguide 
Model (DWM), and some modifications over it. The fields, energy distri- 
bution in all the regions of the structure are obtained. Q-factor has been 
found using the basic definition of Q, based on losses and energy relations. 
Inter-Resonator coupling between two identical DRs for edge coupled and 
broad-side coupled cases have been found, using electric dipole model having 
capacitive coupling. 

Results obtained in this thesis, are in good agreement with the pub- 
lished results, where-ever available. 
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Chapter 1 


Introduction 


1.1 Introduction 

In the present age of miniaturisation of microwave components and 
their applications in millimeteric wave ranges, ‘Dielectric Resonators’ have 
found a prominent place, due to their spectacular properties of small size, 
high Q, temperature stability, high dielectric constant etc. [1], and by now 
these are so well known that it does not need a basic introduction. More 
basic details can be found in [1], [2], [3], [^]. 

Numerous analysis and synthesis methods are existing in various ref- 
erences [2], for characterisation of DRs, like Dielectric Waveguide Model 
(DWM), Mode matching, Integral equations and Green’s function method, 
Finite element method (FEM), Finite difference time domain (FDTD) method 
etc. Almost all of these methods are of rigorous nature and are computation 
extensive. Also, ironically, most of these are dealing with only DRs excited 
in TE 0lS mode, being generally the dominant mode, and thus TMo^s mode 
has been forsaked. 


Off-late, with the increasing popularity of DR operating in TM 07 <j 



mode in' various applications like microwave filter design etc, it’s analysis 
and characterisation has become a challenging task and only a few researchers 
have worked in this field so far. This thesis studies the circuit properties of 
DR operating in TM 0 is mode and presents a simple approximate method, 
which gives reasonably good results in determination of resonance frequency, 
energy distribution, Q- factor and coupling between two adjascent DRs in var- 
ious environments and it’s applications into design of an efficient microwave 
filter. 

1.2 Subscript Notation 


The modes of the dielectric resonators are denoted by HEM mnp 
modes. Here first two subscripts denote the dielectric waveguide modes. For 
m = 0 (i.e. circular symmetry), the fields are split into transverse electric 
to z (TE to z) and transverse magnetic to z ( TM to z) fields, completely 
independent of each other and with no (^-variation. For all m > 0, the fields 
are hybrid, or HEM, having both TE and TM parts. 

The third subscript p denotes an integer number of half wavelength 
variations of the field, as a function of co-ordinate z. The notation TEois 
or TM 0 is, denotes lowest circular symmetric mode of that category (TE or 
TM) (with fields having no variation in <j> direction), having field variations 
of less than one half-wavelength variation within the resonator length (i.e. 
height of a DR disc). 5 therefore signifies a non-integer number, smaller than 



1.3 Organisation of Thesis 


This thesis has been devided into seven chapters. Chapter2 discusses 
the DWM method for calculating resonance frequency, field distributions and 
coupling between two adjascent DRs in TMq u mode. Chapter3 presents the 
electric and magnetic energy distributions in all the regions of the DR’s envi- 
ronment. In chapter4, the Q-factor of DR in various environments, has been 
presented. Chapter5 presents the analysis of DR in a generalised structure 
and determination of inter-resonator coupling. Chapter6 gives the design of a 
filter using DR in TMois mode. Following this chapter comes the conclusion 
and future scope of these studies. 



Chapter 2 


Calculation of Resonance 
Frequency 

2.1 Methods of Analysis 

A number of techniques for analysis of dielectric resonators have been 
reported. Most of these are rigorous analysis methods, which are very 
complicated and computation extensive. Some of these are Mode Match- 
ing Technique, Method of Moments, Finite Element Method (FEM), In- 
tegral Equations and Green’s Function techniques, Finite Difference Time 
Domain (FDTD) etc. Amongst the approximate analysis techniques, Dielec- 
tric Waveguide Model (DWM) is comparatively simple and takes reasonably 
loss computation time. In most of the published literatures on these tech- 
niques, the analysis for only TEqu mode has been discussed, being most 
popular mode for various Microwave applications. However off-late, it has 
been shown that TMou mode of dielectric resonator also shows very inter- 
esting properties. Therefore, now the attention has been diverted to the 
analysis of this mode which was forsaked earlier. In this thesis we have used 
the DWM method for the analysis of T Mqu mode of dielectric resonator. A 



similar analysis for TEqu mode, has already been presented in [6]. 


2.2 Dielectric Waveguide Model (DWM) Method 

In DWM method, a dielectric resonator is considered as the dual of a 
metallic cavity. The DWM method is similar to Itoh’s method [4], except in 
field expressions selected for axial electric field ( E z ). As a truncated hollow 
waveguide becomes a resonant cavity, in the same way the truncated dielec- 
tric rod waveguide becomes a dielectric resonator. For circular symmetric 
configurations, suitable field expressions can be found for axial field compo- 
nents (E. and H z ), from the knowledge of field analysis for dielectric rod 
waveguide. When we apply proper boundary conditions these field expres- 
sions give characteristic equation, the solution of which gives the resonance 
frequency. 

2.2.1 Maxwell’s Equations for TM$ n Case 

For a source free, isotropic and homogeneous medium, the fundamen- 
tal Maxwells equations can be written as:- 

V x E = —jivfiH 

V x H — ju>eE 

V • E = 0 

V • H = 0 (2.1) 

From the equations (2.1) we can easily find the well known vector Helmholtz 
equation :- 


V 2 £ + k 2 E = 0 
V 2 # + k 2 H = 0 


( 2 . 2 ) 



where k is called the wave number of the medium and in general, it is a 
complex quantity. For a dielectric medium of relative dielectric constant e r , 
the wave number k is given by 

k 2 = kle r (2.3) 


where, k 0 = 2tt / 0v /hoeo and / 0 is the resonance frequency. To obtain the 
fields in all the regions of a structure, it is convenient to start from the axial 
component of the electric or magnetic field. Other field components can be 
obtained from this axial field with the help of Maxwell’s equations. For TM 
to 2 mode (be. H z = 0), we can start from the axial electric field be. E z . 
Then vector Helmholtz equation (2.2) reduces to scalar form as follows 


V 2 E z + k 2 E z = 0 


(2.4) 


This equation can be solved by well known separation of variables method. 
The field components for TM to z circular symmetrical case ( d/dj> = 0) 
are expressed as [2] 


E z 

— chosen 

H z 


1 d 2 E z 

H r 

E r 



v 2 dzdr 



= 0 

H* 


0 

0 

jut dE z 
v 2 dr 


(2.5) 


Here, u is defined as the radial wave number, and in general is given by 


z/ 2 = ± ( k 2 — k\) 


( 2 . 6 ) 


where, k z is the axial wave number. 
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2.2.2 Field Distributions in the Resonator 


The cross sectional view of cylindrical resonator placed coaxially in 
a circular metal waveguide cavity is shown in Fig 2.1. For the generality 
of the analysis, the top and bottom circular plates of the circular waveguide 
cavity, can be considered in four possible combinations as follows :- 

1. Top and Bottom circular plates to be PEC 

2. Top and Bottom circular plates to be PMC 

3. Top plate PEC and Bottom plate to be PMC 

4. Top plate PMC and Bottom plate to be PEC 

The resonator material is assumed to be a perfect dielectric charac- 
terised by the real scalar relative permittivity q and the real scalar per- 
meability fi ( = no). To simplify the analysis, this composite structure is 
devided into six regions with a different set of field components. Here the 
proper choice of E z in different regions, depends on the requirements imposed 
by physical geometry and the analyticity of the field at the origin. It also 
should satisfy the wave equation (2.4). Inside the resonator, E z should rep- 
resent standing wave in z direction, and wave number v — k r should be real 
to fulfill the requirement of DR as a energy storing element. The solution of 
equation (2.4) in cylindrical co-ordinates using method of separation of vari- 
ables with above mentioned conditions gives the following form of expressions 
for E- in Region 1 :- 

E z i = {cos((3iz) + B\ sin(Piz)} J Q (k r r ) (2.7) 

where J 0 is the Bessel function of the first kind and zero order, and k r is 
related to axial wave number by separation relation :- 

k 2 r + = e d k 2 (2.8) 
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Outside the resonator E z should be axially symmetric, but decaying in 




PEC/PMC 


PEC/PMC 


y 


PEC circular metallic waveguide 


Figure 2.1: (a) Cylindrical dielectric resonator in a cylindrical waveguide, 
(b) Cross-sectional view of DR in a cavity 

nature in the direction away from the resonator, for confirming the maximum 
energy storage inside the resonator. This decay in radial direction can be 
represented by modified bessel functions, therefore, outside the resonator 
(in the radial direction) the radial wavenumber becomes imaginary. Finally, 
tangential electric fields should vanish at the PEC surfaces and tangential 
magnetic fields should vanish at the PMC surfaces. Again solving equation 
(2.4) with these conditions, E z in the rest of the regions (2, 3, 4, 5, 6) can 



be written as follows. (Note Upper term inside the square braces [ ] is for 
PEC surface at that end of the cavity and lower term is for PMC surface.) 


E z 2 

= b 2 

E z 3 

= b 3 

E z A 

= B 5 

E z 5 

= Be 

Ez 6 

= B r 


cosh{cx a (z - /i 4 )} 

— sinh{a a (z — fi 4 )} 

cosh{a s (z + /i 3 )} 
sinh {a a (z + h 3 )\ 

cosh{o: a (z — /i 4 )} 

— sinh{a a (z — /i 4 )} 

cos(/?i z) + Bisin(; 
cosh{a s (z + h 3 )} 
sinhja^z 4- h 3 )} 


e~ QaflA J Q (k r r ) 
e -a,h 3 j 0 (k r r) 

e -a a h< (/^ 0 ( fcar ) _ BI 0 {k a r)} 
z)} {K 0 (k a r) — BIo(k a r)} 
e-°‘ h * {K 0 (k a r) - BI 0 {k a r)} (2.9) 


where, 


B = K 0 (k a a) / I 0 (k a a) 


( 2 . 10 ) 


I n (x ) and K n (x) are the modified bessel functions of the first and second 
kinds. The separation relations in these regions are 


= kl + a 2 
k 2 T = t s kl + a 2 


k 


2 

a 



k 2 = ~kl + 01 
k\ = -<& ~ a 2 . 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 


The transverse field components in various regions are found by sub- 
stituting (2.7) and (2.9) into (2.5), and are given as 
Region 1 

Eri = {5iC0s(Az) - sin{piz)}Ji{k T r) 

H^i = {cos(Piz) + B\ sin(^iz)} Ji(A: r r) (2.16) 

k r 


9 



be written as follows. (Note Upper term inside the square braces [ ] is for 
PEC surface at that end of the cavity and lower term is for PMC surface.) 


E z 2 

= b 2 

e z3 

= b 3 

E z 4 

= b 5 

e z5 

= B 6 

Ez 6 

= B r 


cosh{o: a (z — h^} 

— sinh{a a (z — h 4 )} 

cosh{a s (z + h 3 )} 
sinh{a 4 (z + h 3 )} 

cosh{a a (z — hi)} 

— sinh{a a (z - h 4 )} 

cos {Piz) + B\ sin(, 
cosh{a s (z + h 3 )} 
sinh{a s (z + h 3 )} 


e~ a ‘ h * J 0 {k r r) 

e -a,h 3 J Q (k r r) 

e ' Qah 4 {K 0 (k a r) - BI 0 {k a r )} 
z)} {I< 0 (k a r) - BI 0 (k a r)} 
e~ a ’ h * { I< 0 (k a r ) - BI 0 {k a r )} (2.9) 


where, 


B = K 0 {k a a)/I 0 (k a a ) (2.10) 

I n (x) and A' n (x) are the modified bessel functions of the first and second 
kinds. The separation relations in these regions are 


= k 2 0 + a 2 a 

(2.11) 

— e s kl + a 2 

(2.12) 

= - <£ 

(2.13) 

= -k 2 o + 01 

(2.14) 

= -e s kl - a 2 s 

(2.15) 


The transverse field components in various regions are found by sub- 
stituting (2.7) and (2.9) into (2.5), and are given as 
Region 1 

Erl = -jr {BiCos(Az) - sin(/3iz)}Ji(/c r r) 
k r 

^ {cosiM + B lS in(Piz)}J x (k r r) (2.16) 

kr 
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Region 2 


E r2 = 

H# 2 = 


Qq-52 -sinh{o: a (z - h 4 )} 

/c r cosh{ct a (z — hi)} 


k r 


cosh{a a (z - h 4 )} 
-sinh{a a (z - /i 4 )} 


e 00/14 Ji(A: r r) 
e -Qa/l4 J x [k r r) 


(2.17) 


Region 3 


Q5-B3 -sinh{a: i (z + h 3 )} 

fc r -cosh{a 5 (z 4- /i 3 )} 


e a,/l3 Ji(k r r) 


H* 3 


ju€oe s Bz 

K 


cosh{a s (z + /i 3 )} 
sinh{a 5 (z + h 3 )} 


e Qi/l3 J\{k r r ) 


Region 4 


(2.18) 


E r 4 = 


a„ -R5 

k'a 


sinh{a a (z — /i 4 )} 
- cosh{a a (z - h 4 )} 


E(j>4 


ja,’e 0 jB 5 

k a 


— cosh{a a (z — h 4 )} 
sinh{o! a (z — /i 4 )} 


e~° th * {K^k^ + BhiKr)} 

(2.19) 


Region 5 

E r 5 = {J5i cos(/?iz) - sin(/? 1 z)}{Rr 1 (A: a r) + 5/!^)} 
k a 

= _7^f£^6 (cos^jz) +R 1 sin(/3 1 z)}{R: 1 (A: a r) + B/ 1 (fc o r)}(2.20) 
k a 


Region 6 


ErG 


a s B 7 sinh{a s (z + h 3 ) } 
k a cosh {a s (z + h 3 )} 


e~ Q ' h * {Ki(k a r) + BI\{k a r)} 


jueg^Eh 

k a 


cosh {a s (z + h 3 )} 
sinh{o; J (z + h 3 )} 


e~*’ h * {Kiikar) + Bh(k a r)} 


(2.21) 


The values of unknown coefficients R ts , can be evalvated by applying the 
boundary conditions at the interfaces of two adjascent regions, and are given 
at Appendix A. 



2.2.3 Resonance Frequency (/ 0 ) 


If k r and Pi are known, then the resonance frequency / 0 can be calcu- 
lated from equation (2.8). The concept of dielectric waveguide model (DWM) 
is used for determination of k r and Pi [7], These are determined as follows. 

2.2.4 Determination of k T 

The resonator is first assumed to extend to infinity in z direction (Fig 
2.2 (a)) for determining k r . We assume that this infinite dielectric cylinder 
carries T Mon mode as the propagating mode, and that k T is radial wave num- 
ber for this cylinder, for r < d. Applying the continuity condition for tangen- 
tial electric field and magnetic field components at the common boundary of 
Region 1 and Region 5, at r = d, gives the following characteristic equation:- 

Jo(A-rd) k a t d Jijkrd) 

{K 0 (k a d)-BI 0 (k a d)} k T {. K x {k a d) + Bh(k a d )} { > 

Radial wavenumber k r is related to k a by :- 

k a = a/(q- 1.0)A:2 _ fc r 2 (2.23) 

k r and k a can be computed by solving equation (2.22) and (2.23). 


2.2.5 Determination of j3\ 

The structure of Fig 2.1(b) is assumed to be a radial double slab 
waveguide (Fig 2.2(b)), alongwhich the TM 0 u guide mode is propagating. 
In other words it is assumed that the axial -wave number Pi is the same 
as that of a double layer slab guide structure, obtained by extending the 
resonator to infinity in the radial direction. When we apply the boundary 
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Figure 2.2: Determination of k r and Pi [7] 


conditions at the interfaces of these regions at 2 : = ±h, 0 < r < d, we get 
one more trancendental equation given by 


777-1 tzm(Pih) — 1.0 1.0 — m 2 tan(/?ih) 

777,! + tan(Pih) m 2 + tan(/?i/i) 


(2.24) 


where, 


??7i 


P 1 


?n 2 


PlCs 

td&s 


coth (a a /i 2 ) 
tanh(a a /i 2 ) 

coth(a s /7i) 

tanh(a,/7i) 


(2.25) 

(2.26) 


We get the value of fc r by solving trancendental equation (2.22) in 
conjunction with equation (2.23). Then we can solve equation (2.24) for res- 
onance frequency /„ for specified height of the DR (i.e. 2h ), after substituting 
the values of p u a«, oc s from the separation relations (2.11) - (2.13). 
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2.2.6 EDC Approach 


As it. will be shown in the next section, the DWM method gives 
slightly higher results for the resonance frequency /o. This error is very large 
when a PEC surface is very close to the circular face of the DR. Mongia 
[7] suggested effective dielectric constant approach (EDC) for calculating the 
radial wavenumber (k r ) for TE 0 is mode, which increases the accuracy of 
DWM method. As per this approach radial wavenumber k T is assumed to 
be the wavenumber- of an infinitely long cylindrical waveguide having same 
radius as that of DR, but having a dielectric constant e e /y, which is less than 
that of the resonator. The relations used in [7] are :- 

ka = {(teff - 1-0) fc 0 2 - fc r 2 } 1/2 (2.27) 

where, 


e eff = (Ujfa + teffb) / 2-0 


£e//a — 


^d — (^d — £effo)hi/d 

e eff 0 


if (h\ < d) 
otherwise 


teffb 
( eff 0 


f Cd ~ (^d — c e //o)h 2 /d 

\ e e//0 


ky(e d )/kl 


if(h 2 < d) 
otherwise 


2.2.7 Numerical Implementation 

Above analysis is very general in nature, and it is applicable for iso- 
lated DR (2a — >• oo, hi -»■ oo,h2 ->• oo), dielectric post having top and 
bottom circular plates PEC (hi « 0, h2 ^ 0), as well as DR in MIC environ- 
ment. Based on this analysis, a menu driven software package (DRTM6) has 
been developed for the calculation of resonance frequency for top and bottom 
circular plates (Fig 2.1) being PEC or PMC. The structural parameters in- 
cluding the dimensions of DR, are assumed to be known. In this program the 



computation is started by finding the root of equation (2.22) using Newton 
Ralphson method with following initial approximate value of k r (e d ){ 7] 


K{(d) = 


0.951x01 + 0.222 {(e d - 1 )(k 0 d) 2 - 0.951xm } /d 


(2.28) 


where, Xoi = 2.405, and is the first root of the equation Jq(x) = 0. Other 
eigen values and the resonance frequency of the resonator can be computed, 
when Ay and k a are determined. Repeated iterative procedure is used to get 
the value of / 0 within tolerable accuracy. 


2.2.8 Numerical Results 

DR Post 

The results obtained for / 0 of a DR post (top and bottom circular faces 
of DR covered with PE(?), are presented in Table 2.1 alongwith the published 
results of rigorous analysis methods of mode matching technique [8] and 
Finite Element Method (FEM) [9]. A comparison of resonance frequencies 
obtained by present method, and that of rigorous analysis method shows 
good agreement. 

Isolated DR 

The results obtained for / 0 of an isolated DR, are presented in Table 
2.2 alongwith the published results and it is found that the present DWM 
method gives very close results on slightly higher side, to those published for 
rigorous methods and EDO method gives very close results on slightly lower 
side of the published results. The error is within 2% and it further reduces 
as value of e d increases. 
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Table 2.1: Comparison of Resonance Frequencies from DWM method, with 
published data for DR post 

hi = h 2 = 1.0 x 10~ 5 mm, = 37.6, e s = 1.0, d = 10.0076mm, a = 12.7mm 


2 h 

(mm) 

Resonance Frequencies (/ 0 ) GHz 

A [8] 

Mode Matching 

/op] 

FEM 

fo 

Present 

13.97 

3.363 

3.38 [15] 

3.389 

11.55 

3.582 

3.613 

3.6109 

10.0076 

3.802 

3.824 

3.8205 

8.95 

4.013 

4.020 

4.0165 

8.17 

4.189 

4.205 

4.2016 

7.56 

4.403 

4.381 

4.3790 

7.076 

4.550 

4.5501 

4.5466 

6.67 

4.710 

4.713 

4.709 


Table 2.2: Comparison of Resonance Frequencies from DWM and EDC meth- 
ods, with published data for isolated DR 


hi = /i 2 = 100mm, e 5 = 1.0, a = 100mm 




d 

Resonance frequencies (/ 0 ) GHz 

fo : 

[Reported] 

fo 

Present (DWM) 

fo 

Present (EDC) 



4.265 

9.44 [10] 

9.603 

9.375 



5.250 

7.5384 [11] 

7.643 

7.415 

-do- 

-do- 

- do - 

7.524/7.6 [12] 

-do- 

-do- 

-do- 

-do- 

- do - 

7.6 [13] 

- do- 

-do- 

1 

5.62 

6.415 

6.1328 [14] 

6.256 

6.069 

79.7 

4.51 

5.145 

5.4072 [14] 

5.419 

5.318 






DR in a Cavity having PMC walls 


The results are obtained for / 0 of a DR placed in a cylindrical cavity, 
having top and bottom circular plates made of PMC, and are shown in Fig 
2.3. These have been compared with the results published in [16] and they 
are found to be in close agreement. The error slowly increases as the PMC 
walls are taken away from the DR. However, the EDC method gives slightly 
low but more closer results. 



/i 2 =15mm, ed= 24, e s =l, d=3. 635mm, a=5.45mm, 2h=4.04mm 


Figure 2.3: Comparison with published data for DR in a cavity with top and 
bottom walls PMC 


DR in a Cavity having PEC walls 

The results are obtained for /o of a DR placed in a cylindrical PEC 
cavity and are shown in Fig 2.4. These have been compared with the results 
published in [16], and it was found that the present DWM method gives very 
high resonance frequency, especially for close distance of PEC circular plates. 



The error reduces as the PEC walls are taken away from the DR surface. The 
EDC method also somewhat tends to fail for close proximity of PEC walls 
and thereafter it gives a feel of fall in resonance frequency with increase in 
wall distance, till it saturates for an isolated DR case. The reason for such a 



Substrate Thickness (hi) (mm) 


/?. 2 =15mm, e d =24, e s =l, d=3. 635mm, a=5.45mm, 2h=4.04mm 

Figure 2.4: Comparison of results from pure DWM and EDC methods, with 
published data for DR in a PEC cavity. 

high results obtained from DWM method are analysed to be happening due 
to mismatching of the fields in the outer regions of the cavity e.g. Region 
2-4, 3-6, 4-5, 5-6. Therefore for PEC cavity we had tried to use the modified 
DWM method. 


2.3 Modified DWM Method 

If we match the boudary conditions for continuity of tangential fields 
at, r = d, h < z < h A and -h < z < -h 3 in Fig 2.1 then we find a new 




trancendental equation as 


Jo(k r d) k a .h(k T d) 

{Ko(kad) - BI 0 {k a d)} k r {K\(k a d) + BI x (k a d)} 


(2.29) 


Similarly, on matching the tangential fields at z = ±h, d < r < a we again 
find a new trancendental equation as follows 


tan(/3j/i) — 1.0 1.0 — m' 2 tan(/?i/i) _ 

m' 1 +tan(^ 1 fi) m 2 + tan(j3 x h) 

where, 


(2.30) 


m 2 


Pl_ 

®a 

0l£s 

a s 


coth(a a /i2) 

tanh(or a /i 2 ) 

coth(a s /ii) 

tanh(a s /ii) 


(2.31) 

(2.32) 


In order to include the effects of these additional trancendental equa- 
tions (due to fields in the outer regions) we modify the trancendental equation 
to be solved to determine the resonance frequency, by combining equation 
(2.24) and (2.30), as follows 


/ 77? i tan ((3 x h) - 1 1 - m 2 tan(/? I /t) 

\ 7??,i + tan(/?x h) tt7 2 + tan((3 x h) 

/ m[ tan (^i ft) - 1 _ 1 - tan(/? x fe) \ _ Q ^ ^ 
\ 777i + tan(/?i h) tt 7 2 + tan(/?ifi) J 

Therefore, for a particular position of PEC walls either first term is 
zero, or second term is zero, or both are individually zero. The effect of 
equation (2.29) has however been neglected, as any change in this equation, 
changes all the results drastically. 



2.3.1 Numerical Implementation 

Now, solving equation (2.22) in conjunction with equation (2.23) for 
radial wavenumber k r , and equation (2.33) in conjunction with equation 
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(2.25), (2.26), (2.31) and (2.32), we can find the resonance frequency for 
a DR placed in a PEC cavity, as done in section 2.2. 

2.3.2 Results for DR in a PEC Cavity 

The results are obtained for /o of a DR placed in a cylindrical PEC 
cavity, using modified DWM approach and are compared with that in [16] 
and it was found that the new resonance frequency is in close agreement, 
especially for very close distances of PEC from DR’s flat surface, however 
the error increases slightly as the PEC wall is taken away from the DR, as 
shown in Fig 2.5. 



/i2=15mm, ed= 24, e,=l, d=3. 636mm, a=5.45mm, 2h=4.04mm 


Figure 2.5: Comparison of results from different theories, with published 
data for DR in PEC cavity 
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2.4 Final Selection of Theory 


On perusal of the results obtained from various theories used, it was 
decided to use Pure DWM or Pure EDO for the case of DR in a cavity having 
top and bottom walls PMC. However, for DR in PEC cavity, upto a distance 
of 0 < h\/d < 0.2 and/or 0 < h^d < 0.2, modified DWM is used and 
thereafter we switchover to the modified EDC i.e. using e e//0 for hy/d > 0.2 
and/or h^/d > 0.2, for calculation of radial wavenumber k T . The results of 
this, are shown in Figure 2.6, in which the sudden jump in the variation is 
at the place where theory used is switched from Modified DWM to Modified 
EDC. 



/? z =15mm, u= 24, e,=l, d=3.635mm, a=5.45mm, 2h=4.04mm 
Figure 2.6: Variation of Resonance frequency with hy using modified method. 
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2.5 Effect of Design Parameters on /o 


Various electrical and physical parameters of the circuit, which are 
under control of the designer at the time of fabrication of the circuit, may 
change the resonance frequency of the DR in the cavity. Informations pre- 
sented in the Figs 2.6 to 2.10, shows the impact of these control variables 
over resonance frequency / 0 . 



/i! =6. 98466mm, /i 2 =12.318394mm, e,=l, d=8.88957mm, 
a=12. 69937mm, h=6. 0957mm 

Figure 2.7: Variation of Resonance frequency with dielectric constant of DR 

M 


2.5.1 DR in a PEC cavity 

The effects of separation (h 2 ) between top PEC plate and the top 
surface of the DR, over f 0 , are similar as that of h x as shown in Fig 2.6. 
Variations of resonance frequency with h 2 can be explained by the cavity 
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pert.orbation theory i.e. when a metal wall of a resonant cavity is moved 
inwards, the resonance frequency of the cavity will incrcafio if the stored 
energy of the displaced fields, is predominantly electric, which is the case 
with TMois mode. 



/?i=1.0xl0 -5 , /i 2 =20.319mm, e s =l, d=6. 34969mm, a=12.69937mm, 

2h=5. 07975mm 

Figure 2.8: Variation of Resonance frequency with dielectric constant of DR 

(<*) 


Fig 2.7 shows the effects of the dielectric constant (e d ) of the DR, 
on it’s resonance frequency. We find that the resonance frequency decreases 
with the increase in the value of e d . Hence the material of the DR should be 
temperature stable to get stable resonance frequency. Here, modified EDC 
part of the theory has come into play and gives close results to those reported 
in [24]. Fig 2.8 also shows the variation of e d with / 0 but for a DR placed on 
the PEC surface. Here, modified DWM part of theory has come into play and 
gives very close results to those in [24]. Fig 2.9 shows the effects of substrate 
dielectric constant (e s ) over the resonance frequency. It is found that the 
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Substrate Dielectric Constant (e s ) 
h 2 = hi=3mm, d=3. 635mm, a=5.45mm, 2h=4.04mm 


Figure 2.9: Variation of Resonance frequency with dielectric constant of sub- 
strate (e 5 ) 

resonance frequency decreases with the increase in the substrate dielectric 
constant. Fig 2.10 shows the effects of height of the DR (2 h) over / 0 , since 
the DR height can also be used as a prime design parameter in excitation of 
a particular mode of the DR and achieving a particular resonance frequency. 
With increase in DR height the resonance frequency is found to be decreasing 
rapidly. 


2.5.2 DR in a Cavity with PMC walls 

The effects of separation (/i 2 ) between top PMC plate and the top 
surface of the DR, over resonance frequency / 0 , is found to be decreasing 
with the increase in the separation /i 2 , and is same as that of hi as shown in 
Fig 2.3. The variations of resonance frequency with the dielectric constant 
of the DR (e rf ), substrate dielectric constant (e*), DR height (2 h) are found 
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Figure 2.10: Variation of Resonance frequency with height of the DR (2 h) 


to be similar to the case of PEC cavity. 


2.6 Coupling Between Two DRs 


For the TMqu mode, coupled identical resonators placed in the con- 
figuration shown in Fig 2.11, inter-resonator coupling can be computed using 
the resonance frequencies for a DR placed in a PEC cavity (f S h) an d that of 
placed in a cavity having top and bottom circular plates to be PMC (/op). 
Kobayashi etal. [16] has reported this method of computing the coupling co- 
efficient. The calculated results show that fop > f sh ; which is in contrast with 
the usual inductively coupled case for the TEou mode in which f„h > fop- 
Colm [3] had shown that the inter-resonator coupling is caused by an axial 
magnetic dipole for the TE m mode. By analogy with this case, we expect 
that, it is done by an axial electric dipole for the T M 0 u mode. Then we shall 
consider the capacitively coupled LC resonant circuit shown in Fig 2.12(a). 
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Fig 2.12(b) shows it’s equivalent circuit, where f op and f s h are the resonance 
frequencies when the symmetric plane T(T') is open circuited (PMC) and 
short circuited (PEC), respectively. From the Fig 2.12(a) and (b), we have 




Figure 2.11: (a) Circuit configuration for inter-resonator coupling, (b) DRs 
having air in-between, (c) DRs having substrate in-between 


/o 


I 

2ttVLC 


(2.34) 


where / 0 is the resonance frequency of the LC resonant circuit (from Fig 
2.12(a)). 


fop /o 


(2.35) 
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Figure 2.12: (a) Capacitively coupled LC resonant circuit, (b) Equivalent 
circuit. 


fsh 


Jo 



(2.36) 


where k is the inter-resonator coupling coefficient and is defined as [16] 


C + C m 

From equations (2.35) and (2.36), it can be shown that 
f 2 - f 2 

j J op J sh 

~ 1^+JTh 


(2.37) 


(2.38) 


2.6.1 Numerical Implementation 

f and f, h ran be found for both the structures shown in Fig 2.11(b) 
& (c), and therefore the inter-resonator coupling coefficient can be computed 
using equation (2.38). Fig 2.13 shows the comparison of computed results 
as against published data [16] for the structure of Fig 2.11(b) without any 
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substrate (i.e. e s — 1.0), and neglecting the effect of dielectric support (e s = 
1.03) used in [16]. Here it can be observed that the present theory gives 



hi = h.2=3mm, 6^=24, e s =5, d=3. 635mm, a=5.45mm, 
2h=4.04mm 

Figure 2.13: Comparison of computed inter-resonator coupling with pub- 
lished data in [16]. 

slightly higher results for inter-resonator coupling coefficient for very small 
air gap and very low results thereafter, which is due to the slight error in the 
computation of f op (being slightly high by pure DWM method) and f,h (being 
slightly low by the modified DWM method). However, the result achieved 
is that, the inter-resonator coupling reduces rapidly as the distance between 
the two dielectric resonators increases, and that they couple electrically with 
each other when the dominant evanescent waveguide mode is the TM mode. 




Chapter 3 


Energy Distribution 


The relative energy distribution in different regions of the circuit 
structure involving dielectric resonator, is very useful in deciding how to do 
the coupling of DR with the other components. This chapter will deal with 
the electric and magnetic energy distribution for TM 0i s mode at resonance. 


3.1 Electrical Energy Distribution 

Total electrical energy stored in the structure given in Fig 2.1, is the 
sum of individual electrical energy stored in each region. It may be calculated 
from the field expressions derived in chapter 2. 

it; = E "» < 31 > 

i=i 

where, IT,., is the electrical energy stored in the i th region of the structure 
occupying volume V and having dielectric constant e r j. W e i is given by 

IF,,- = Ei-E‘ dv (3.2) 

4 Jv 
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Wlioro, 1+1 indicates the complex conjugate. For TM m mode only E z and 

E r components of electric field exist, therefore for different regions we get 
stored electrical energy as follows 
Region 1 

^0 fd r27r th . ^ 

U ' el = ~T Jo J 0 J_ h { E *-E'zi + 4i*4i) rdrd<j>dz (3.3) 
Region 2 

]le2 ^ J Jo J 0 J h (■ E * 2 ' E z2 + Er2 • E; 2 ) r dr d<j> dz (3.4) 
Region 3 


^ 7f, 3 — ^ Jo Jo J_ hz i Ez3 ' E z3 + E r3 ' E* 3 ) v dr dj> dz (3.5) 
Region 4 

W ,r e4 = J d J 0 J h (Ez 4 'E* z4 + E ri -E*^j rdrd(j>dz (3.6) 

Region 5 

ir e5 = | JJ Jj* fj (4 5 • e ; 5 + E r5 ■ e; 5 ) r dr dj> dz ( 3 . 7 ) 

Region 6 

n; 6 = ~ JJ jj* jj (4 6 • 4*6 + Ere • 4e) r dr d<f> dz (3.8) 

After substituting the values of fields from chapter 2, and solving integrals, 
we get the expressions for these electrical energies as follows 
Note Here the upper term in square bracket [ ] corresponds to correspond- 
ing circular plate of cavity to be PEC and lower term corresponds to the 
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corresponding plate of the cavity to be PMC. 


IP, 


el 


n; 2 


w. 


e3 


h; 4 


Ues 

ir e6 


€d | ^16 Ijo + 01 -^16 -4x234 j 


A 0 Bf | 

^(h2,/l4,C*a) 

. ^2,^4, <*a) 

h(3 + C*a 

Ao B%e s | 

ft 

A^l>h3,Q 5 ) 

ho + o? 3 

A 0 Bl | 

^(/l 2 ,/l4,Cta) 

. ^2,^4, a a ) 

7/c/o + ocl 

A 0 Bj { A 

L 16 IkIO + Pi A' 16 A 5678 } 

AqBjCs | 

. t(hi,hz , q s ) 

IkIO + a 2 s 


^(/i2,/i4 ) aa) 

t* 

(h2Ai,a a ) . 
ft 

l (hi i hz,a s ) m 
Hh2,h4 ,Q a ) 

. ^(/l2,^4, Cta) 



Hhuhz,a a ) 

ft 

. L {hi,hz,a a ) 



(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


where, 


Ao = 
A 1234 

^5678 


COTT 


[ej_ 

k 2 r 

Ieki 

*2 


(3.15) 

(3.16) 

(3.17) 


Here constants B is are same as in chapter 2, and are given at Appendix A. 
Other constants are given at Appendix B. 


3.1.1 Electric Energy Filling Factor 

Energy filling factor is an important quantity used in connection with 
the analysis of partially filled cavities. This is defined as the ratio of the elec- 
tric energy for a given region to the total electric energy stored in the system. 
Based on the above analysis for electric energy, a program has been devel- 
oped to find the energy in all the regions of the structure. Results obtained 
are tabulated in terms of energy filling factor for each region for different 
dielectric constants of the resonator and the substrate. It can be observed 
from the Table 3.1, that for e d = 38.0, more electrical energy gets confined 
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Table 3.1: Electric energy filling factor for various regions of resonant struc- 
ture 

d = 3.63577i7n, a = 5.45mm, /t t = 2.0 mm, = 2.0mm, 2 h — 4.04mm 



e s 

Energy Filling Factor (= W ei /W e )% for region No. 

1 

2 

3 

4 

5 

6 

24.0 

2.22 

51.96 

0.476 

1.156 

6.618 

24.58 

15.20 

30.0 

2.22 

54.50 

0.383 

0.914 

6.225 

23.76 

14.20 

30.0 

10.0 

34.12 

0.209 

3.933 

3.870 

13.45 

44.40 

38.0 

2.22 

56.79 

0.305 

0.715 

5.879 

22.96 

13.34 


within the resonator than when e d =24.0, therefore energy leakage in other 
regions increases when dielectric constant of DR (ej) decreases. Similarly, 
using higher dielectric constant substrate, energy contents in the region 3 
and 6 increases. Also, since about 40-50% electrical energy is distributed 
in outer regions, electric coupling of the DR (in TM 0 is mode); with other 
circuit elements is very strong. 


3.2 Magnetic Energy Distribution 

We can also find the magnetic energy stored in the structure in Fig 
2.1, in the same way as we have done for the electric energy. The net magnetic 
energy is the sum of magnetic energies of individual regions. 

It™ = E U'„, (3-18) 

t=l 

where, IT,,,; is the magnetic energy stored in the i t h region of the structure 
occupying volume V and having relative permeability /x ri . W mi is given by :- 

H„i — f Hi ■ dv (3-19) 

4 Jv 
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For T Alois mode only H# component of magnetic field exists, therefore for 
different regions we get stored magnetic energy as follows 
Region 1 


Li o f d fh , \ 

^ ml = -J J Q J o J_ h ( H# ■ ) r dr d<j> dz (3.20) 

Region 2 

Un fd- r2ir rht . _ _ . 

11 m 2 = 1 J o J o j h ( H <t>2 ■ H$2 ) rdrd(j)dz (3.21) 

Region 3 

n " m3 = T Jo C J-h (^ 3 • ^3 ) r dr d( t> dz (3.22) 

Region 4 

W n4 = ^ J* fj (H^ -H^) r dr d<j> dz (3.23) 

Region 5 

w mi = ^ jT J** f h (H # 5 - H; 5 ) r dr d<f> dz (3.24) 

Region 6 

in* = j f‘ [' J2 {&» ■ HU) r dr d<fr dz (3.25) 


After substituting the values of fields from chapter 2, and solving integrals, 
we get the expressions for these magnetic energies as follows 
Note Here the upper term in square bracket [ ] corresponds to correspond- 
ing circular plate of cavity to be PEC and lower term corresponds to the 
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corresponding plate of the cavity to be PMC. 


HVn, 

ir m2 


A )^d { ^16 j 4 j 234 } 


11 


m3 


D 0 B% 


DoBje 2 s 


(/l2,/l4>Ct a ) 

L HhtJi 4,a a ) 

f 

L {huh3 y a 3 ) 
(hl,/l 3 ,Q 5 ) 


t 


w mA 

Urns 

^ m6 


A>sf 


ft 

(/i2,/i4,a a 

L Hhito&a) J 

DqB\ { j4i6 ^5678 } 

t'< 


DoBhl 


{hi,h3,ct,) 

L kh iMv.) J 


^1234 

^■1234 

^5678 

•^5678 


where, 


D n 


A 


1234 


■^5678 = 


Mqtt (w 0 eo)" 


Iej 

k r 

Ieki 

k2 a 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 


(3.32) 

(3.33) 

(3.34) 


Here constants B{ s are same as in chapter 2, and are given at Appendix A. 
Other constants are given at Appendix B. 


3.2.1 Magnetic Energy Filling Factor 

This factor for a region represents the part of total magnetic energy in 
that region. For TM 0 is mode, though it is not of much importance because 
most of the magnetic energy is contained within the dielectric resonator and 
almost negligible amount is in the other regions, therefore magnetic coupling 
is extremely weak for TM 0 is mode. Table 3.2 shows the magnetic energy 
distribution in all the regions of the structure. We can observe from Table 3.2 
that as was expected most of the magnetic energy (90-95%) is stored within 
the DR and only fractional amount (5%) is distributed in the other regions. 
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Table 3.2: Magnetic energy filling factor for various regions of resonant struc- 
ture 

d = 3.63577177?, a = 5.45mm, hi = 3.0mm, = 3.0mm, 2 h = 4.04mm 




Energy Filling Factor ( 

= W mi /w m )% for region No. 

1 

2 

3 

4 

5 

6 

24.0 

2.22 

91.12 

0.029 

0.166 

0.791 

4.289 

3.605 

30.0 

2.22 

93.354 

0.018 

0.100 

0.576 

3.283 

2.667 

30.0 

10.0 

70.793 

0.010 

2.780 

0.378 

2.840 

23.19 

38.0 

2.22 

95.070 

0.011 

0.059 

0.417 

2.482 

1.957 


As the dielectric constant of the DR increases the stored energy within the 
DR increases, therefore, magnetic energy does not play important role in the 
coupling, when operating in the TMqu mode. 
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Chapter 4 


Unloaded Quality Factor (Q 0 ) 


4.1 Quality Factor 


Quality factor of a circuit is a very important parameter because it 
limits the overall performance of the circuit. Higher value of Q is required in 
design of narrowband filters. The unloaded Quality factor of a DR placed in 
an MIC environment or a PEC cavity, depends upon losses in the DR’s mate- 
rial, conductor loss in the shielding walls and dielectric loss in the substrate. 
Many rigorous analysis methods have been reported to find the Q-factor 
which are quite accurate, but are highly complicated and numerical exten- 
sive. Frequency perterbation technique has been reported for computing the 
Q-factor for TE 0 \s mode by Kajfez [17]. A similar approach can however be 
derived for computng the Q-factor for TMou mode, for finding Q-factors due 
to dielectric rod (DR) loss, dielectric support loss, but it cannot be used for 
finding the Q-factor due to conductor losses [16]. Although this method is 
quite accurate, lot of caution may be required while applying this method, 
because the relative change in resonance frequency due to slight perterbation, 
is very large for close PEC walls and causes error. 

In this thesis Q-factor is determined for the structure of Fig 2.1 by 



computing the stored energy and power loss. This method is approximate 
but, is quite accurate and data calculated by this method can be used for 
design purposes. 


4.1.1 Definition of Q-factor 


The Q-factor is assurance of the performance or quality of resonator, 
and is a measure of energy loss or dissipation loss per cycle as compared to 
the stored energy. Q-factor is defined by the well known relation :- 

„ . max energy stored 

Q = 2 7T -T- - - , r (4.1 

avg. energy dissipated per cycle 

Let U'e and W m are the stored electric and magnetic energies respec- 
tively, /o is the resonance frequency and P is the total loss, then Q will be 
given by :- 


Q = 


(H'« + Wrn) 

P/I 0 


(4.2) 


4.2 Determination of Conductor Losses 

The tangential magnetic field generates the conductor current on the 
top and bottom metal plates and cylindrical side plates of the waveguide 
section. Let J c is the surface current density on the metal surfaces, then the 
conductor loss P c can be expressed as :- 

P c = l R s J\l \ 2 ds (4.3) 

where, R s is the real part of surface impedance and is given by 


R s = \Af /o Mo / o 


(4.4) 



where, o is the conductivity of the metal surfaces, and J c = z x H, and z 
is the normal vector to the surface. 


4.2.1 Losses in Top Circular Plate 

The conductor losses in the top circular plate of the cavity can be 
given as 

Pc.{top) = 2 R s J I Jc(top) I ds 

In TMq\s mode only H $ component of the magnetic field exists, hence 
Jc(top) = Z x H# (for region 2 & 4). 

Region 2 

J c 2 {top) = z x. Hfpi (4.5) 

e- 2a ' h ‘Cf‘ Jl(krr)vdrd4, 

=>■ P C 2{top) = nR* (weo B2) 2 e 2ctah *Ai2M ( 4 - 6 ) 

Region 4 

del (top) = z x Hfi 4 (4T) 

Pci(top) = 6-200,14 Jo Id [KiM + B h{k a r)f r dr d<i) 

=£■ Pci(top) = 7T R s (lucqBs) 2 e 2aahi A^Qjg (4.8) 

4.2.2 Losses in Bottom Circular Plate 

Region 3 

= \ R ,^i^s:i: Jim r^i * 

=s. P mM) = nR,(u«xM 2 e- 2 ° M A I23 , (4.9) 
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Region 6 


■Pc6(6oi) — 2 ^ s ( ) e 2a ’ h3 J 0 J d [K\{k a r) + BIi(k a r)] 2 r dr d(f) 

=> Pc6(bot) = kR* (we 0 e 3 B 7 ) 2 e~ 2a ‘ h3 A 567S (4.10) 

4.2.3 Losses in Cylindrical Side Walls 

The conductor losses in cylindrical side walls of the cavity, will be due 
to Hj component of magnetic field in Region 4 , 5 and 6, and J c(3i£fe) = rxH$. 

Region 4 

P cA(side) = 5 j e~ 2 Qah *KImJ^ [cosh {a a {z - h 4 )}] 2 dr dz 

^ Pci{side) = TT a R s j Kim t'(h 2 ,h 4 ,aa) (4-11) 

Region 5 

Pc5(side) = ~Rs (~~ ^ 6 ) Kl 456 J 0 j_ h [cos{Piz) + Bismdhz)} 2 dr dz 
=> Reside) = * aR ,(^^ 2 I<ImA 16 ( 4 . 12 ) 

Region 6 

PcHsUe) = ( ‘‘'"^ ~) e~ 2 °‘ l "KI<x j 0 (cosh {a, (z + h 3 ))] 2 a ir dz 

=» P^,uU) = r a R, (^Al}KI tx ( 4 . 13 ) 


where, 

A J456 = [A 1 ( k a ci ) + BI\ (fc a a)] 

B = K 0 {k a a) /h{k a a) 

B is are given at Appendix A and other constants are given at Appendix B. 
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4.3 Conductor Quality Factor ( Q c ) 

If total coductor loss is given by P c then 


P c — (P c 2(<op) + Pc4(top)) + (P =3(feoi) + Pc6(6oJ)) + (P C 4(stde) + Pc5(side) + Pc6(tide) 
=> Pc = P c[top) + P c (bot) + Pc(side) (4-14) 


where, P c (( 0p ), P c (bot), P c (stde) are the total conductor losses in top, bottom 
and side walls respectively. Q-factor due to conductor losses ( Q c ) can be 
defined as 


Qc = 2 7T /o 


(We + w m ) 

" Pc 


(4-15) 


therefore, from equation (4.14) we get, 


Pc Pc(top) . Pc(bot) 

2tt/o (W e + W m ) ~ 2tt/o (W e + W m ) + 2tt/ 0 {W e + W m ) 

1111 
— — -f. -f 

Qc Qc(top) Qc(bot) Qc(side) 


Pc(side) 

27r/ 0 (W e + W rn ) 
(4.16) 


where, Qc(to P ), Q c (bot), Qc(side) are the Q-factors due to conductor losses in 
top, bottom and side walls respectively. 


4.4 Dielectric Quality Factor ( Qd ) 

Intrinsic Q-factor of the DR (Q dr ) is dependant upon the loss tangent 
of the DR material. Q d of DR in MIC environment or in a cavity, does 
not exactly equals l/tan<5. It also depends upon electric energy distribution 
inside and outside the resonator [18]. A similar argument is applicable to 
Q-factor due to dielectric loss in the substrate. By definition Q d is given by 



where, W = \V e + fT m « 2W e , and P d = 2ujW d t&n6. Here, W d is the 
electrical energy stored in the dielectric region [18]. 

Hence, Q-factor due to dielectric loss in different regions are 
Q-factor due to dielectric losses in the DR is 


Qdr 


n; 


11 ci tan S dr 


(4.18) 


Q-factor due to dielectric losses in the substrate (Region 3 and 6) is:- 


Qds 


h; 


(ir e3 + W e6 ) tan S ds 


(4.19) 


Here, tan 6 dr and tan 8 ds are the loss tangents of the resonator and substrate 
respectively. Total dielectric loss in all the regions, is equal to:- 


Pdt — Pdr + Pds 

=*> P dt = 2u >\V e i tan 5 dr + 2oj (W e3 + W e6 ) tan 6 ds 

Pdt _ H'ei tan 5 dr (11^3 + W e5 ) tan 5 d , 

^ 2u AV e ~ \V e + W e 

J_ _ JL JL 

Qdt Qdr Qds 

Wg 

^ dt VHei tan 6 dr + (lT e 3 + W t $) tan 6 ds 


(4.20) 


4.5 Overall Unloaded Q-factor (Qo) 


Once the conductor Q-factor and the Q-factor due to dielectric losses 
are calculated, overall unloaded Q-factor can be calculated as 

_L _ J_ J_ 

Qo ~ Qc + Qdt 


AO 


(4.21) 



Table 4.1: Comparison of Q dr , Q c > and. Q 0 with published data [16] 
d. = 3.635 mm, a = 5.45mm, hi = h 2 = 14.54mm, 2 h = 4.04mm 
(s — 1-0, e d = 24, tan^r = 4x 10~ 5 , a = 52.2 x 10 6 g/m 



Qdr 

Qc 

Q o 

Reported [16] 
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Note Here effect due to support dielectric (e r = 1.031) in [16] has been 

neglected. 

4.6 Numerical Implementation 

Based on the above analysis, a program has been developed to deter- 
mine the overall unloaded Q-factor of a DR in an MIC environment or kept 
in a PEC cavity. The results have been generated to obtain the conductor 
Q-factor (Q c ), dielectric Q-factor (Q dt ) and overall unloaded Q-factor (Q 0 ) 
and these have been compared with published data [16] in Table 4.1. 

Data for Q 0 have also been generated for DR in MIC environment, as a func- 
tion of substrate thickness (hi), with different parameters such as top plate 
distance (h 2 ) , e d , e s etc. and they are shown from Fig 4.1 to 4.5. 

Fig 4.1 shows the variation of Q c with substrate thickness for different 
values of h 2 . It can be seen that, as the conductor plates are taken away from 
the DR, the conductor Q-factor increases, which is due to the decrease in the 
conductor losses in those plates due to decaying magnetic field with distance. 
This effect on Q c causes the increase in the overall unloaded Q-factor, which 
can be seen in Fig 4.2. Fig 4.3 shows the effect of substrate dielectric constant 
i.e. the Q 0 decreases with increase in the substrate dielectric constant due 
to increase in the dielectric losses in the substrate. Fig 4.4 shows the effect 
or DR. dielectric constant over Qo, and that also follows the same analogy 
of decrease in Q„ with increase in the u due to increase in the dielectric 
losses with increase in £ „. Fig 4.5 shows that DR height also plays a role in 
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Figure 4.1: Variation of Q c with conductor plate distance 


controlling the Qo, as with increase in DR height, the resonance frequency 
decreases and thereby decreasing the losses. Therefore Qo improves with the 
increase in the DR height. 
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Figure 4.2: Variation of Q 0 with conductor plate distance 



e rf =24, /? 2 =2.0, d=3.635, a=5.45, 2h=4.04, tan 5ds = tan idr = 4 x 10 5 

Figure 4.3: Variation of Q 0 with conductor plate distance, for different values 
of substrate dielectric constant (e 5 ) 
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e s =2.22, e d =24.0, lx 2 = 2.0, d=3.635, a=5.45, tan** = tan f(ir = 4 x 10“ 5 


Figure 4.5: Variation of Q 0 with conductor plate distance, for different values 
of DR Height (2/i) 





Chapter 5 


DR in a Generalised Structure 


Inter-resonator coupling between two DRs, found in chapter2, gives 
very high (for close PEC) and very low (for far PEC) results. The main 
reason for this error, was the inaccuracy in determination of f sh (resonance 
frequency of DR in a PEC cavity) and thereby the error was incorporated in 
determination of inter-resonator coupling using Kobayashi’s method [16] and 
equation (2.38). In order to determine the inter-resonator coupling for the 
structures shown in Fig 2.11 (b)&(c), using Skalicky’s method [19] and that 
reported by Byzery [20], a much more generalised cavity structure containing 
an axially symmetrically placed DR (Fig 5.1) was concieved. 

The other major advantage which was achieved simultaneously, was 
that the structure of Fig 5.1 could be used for determination of resonance 
frequency of dielectric resonator operating in TM 0 u mode, in practically all 
the possible types of environments like Isolated DR, DR post, DR in MIC or 
suspended substrate environments etc. 



5.1 Resonance Frequency 


The resonance frequency of the structure in Fig 5.1, is found using 
pure DWM method, by deriving the trancendental equations exactly on the 
similar lines as done in chapter 2, for a six region structure. For ease of 
analysis, the structure of Fig 5.1(a) can be devided into ten regions as shown 
in Fig 5.1(b). 

5.1.1 Field Distributions 

Proceeding on exactly similar lines of section 2.2.2, the field equations 
and corresponding separation equations in all the ten regions, can easily be 
written as 
Region 1 

E z i = {cos(/?i z) + BiSin(Piz)} J 0 (k r r ) 

Eri = ” {Bi cos(ftz) - sin{Piz)} Ji{k r r) 

K r 

H t , = - -p— {cos(ftz) + S t sin(i'iiz)} Ji(k r r) 

ft = e d ft - ft (5.1) 

Region 2 

E z 2 = B 2 {cos{(3 2 z) + B 2 sm{p 2 z)} J Q (k r r) 

E r2 = -j^B 2 {Bz cos(P 2 z) - sin {p 2 z)} Ji{k r r) 

fcp 

H t 2 = — tr— - - B2 (cos(^2 2 ) + St sin(/Sz)} Ji(k r r) 

ft = - ft < 5 - 2 ) 

Region 3 

E , = a \ “sh{ Q . 3 (z - he) } 1 UKt) 

23 4 [ — sinh{a 5 3(^ — he)} 



circular waveguide 



Figure 5.1: (a) Cylindrical dielectric resonator in a partially filled cylindrical 
metal waveguide, (b) Cross-sectional view of one stage. 
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cosh{a s2 (z + /i 8 )} 
sinh{a s2 (z + h 3 )} 
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- coshja^z + h 3 )} 
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J 0 {k T r) 


e a,2hs Ji{k r r) 

e -a. 2 h B jj(^ rr ) 
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s2 
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E z 6 

E^ 


= B« 


H(j>e 


a S 3 ry 

-J — -£>8 

k a 

jue 0 e s3 


k a 
-e S 3^o 


cosh{a i3 (z - h 6 )} 
-sinh{o; a 3 (z - h 6 )} 

sinh{a a3 (z - /i 6 )} 
cosh{a j3 (z - /i 6 )} 

— cosh{a s3 (z - he)} 
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Be 


a 


s 3 


e _Q ‘ 3 ' 16 {if 0 (A: ar ) - BI 0 {k a r)} 
e ~ a ‘ 3h 6 {^(W + B/^V)} 

e~ a,3he {K^r) + BIxikar)} 
(5.6) 


Region 7 


£ 27 = Bg {cos (/? 2 z) + B 3 sin(/? 2 z)} {i^o^aO - BI 0 (k a r)} 


Erl = t-B 9 {B 3 cos(/3 2 z) - sin (&z)} {K x (k a r) + BI x (k a r)} 

H *7 = ~—^——-B 9 {cos (/? 2 z) + B 3 sin(/? 2 z)} {K x {k a r) + BI x {k a r)} 

k 2 a = —tsikl + Pi (5.7) 


Region 8 


E zS = B xo {cos {P x z) + B x sm{p x z)} {K 0 {k a r) - BI 0 (k a r)} 

E r8 = ^-B x o {B x cos(P x z) - sin^z)} {K x (k a r) + BI x (k a r)} 

H 4,8 = - ~~ —B x0 {cos{p x z) + B x sin(p x z)} {K x (k a r) + BI x {k a r)} 
K a 

kl = -e ss k% + Pl (5.8) 


Region 9 


E z 9 = B n {cos(/? 3 z) + B 6 sin(/? 3 z)} {K 0 (k a r) - BI 0 (k a r)} 

E r9 = —Bn {B 6 cos(// 3 z) - sin(/3 3 z)} {K x {k a r) + BI x (k a r)} 
k>a 

H„ = {cos(ftz) + B 6 sm(ftz)} {K,(Kr) + Bh(Kr)} 

l'a — + Pi (5-9) 
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cosh{a i2 (z + M} 
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cosh{a s2 (z + /i 8 )} 
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B 12 


a? 


e -a.^ {Ko{kar ) _ J3/ 0 (fc a r)} 

e~ a ‘ ihi {K x (k a r) + BI x (k a r)} 

e -a, 2 hs {^(^O+B/^fcar)} 
(5.10) 


Here, 


B = K,(Ka)/Io(koa) 


(5.11) 



Here also, upper term in the square brackets [ ] refers to the corresponding 
circular plate to be PEC and lower term refers to the corresponding circular 
plate to be PMC. Now, on applying the boundary conditions and equating 
the tangential field components at the boundary of each region and thereafter 
equating the expressions of like field constants, we get 

Jo(M) htd MKd) » ,_ 19 , 

{A'o(M) - BI 0 (k a d)} ^ k r e ss {Ki{k a d) + Bh(k a d)} ~ [ * AZ) 


tdlh. + e, 4 /Am 3 tan(/?i k) _ £d/?3 ~ tan (fak) 

£*4Am 3 - e d (3 2 tan(/?ih) + e d p 3 tan(/?i h) 


where, 


1 + B 3 tan (/?2 A-) 
Bz - tan (/ 3 2 h ) 


(5.14) 


1 - Be tan (/? 3 /i) 
B 6 4- tan (P 3 h) 


(5.15) 


and, 


p 2 ts 3 [ coth(o, 3 /r 4 ) 
e, 4 a, 3 tanh(o, 3 /i 4 ) 


1 + B 3 tan (P 2 h 3 ) 
tan (fohs) - B z 


= mi (say) 


(5.16) 


Pz^s 2 [ coth ( 0 , 2 / 13 ) 

e s iOc s2 I tanh(a s 2 /i3) 


1 — Bq tan (^ 3 / 17 ) 
tan (P 3 h 7 ) + B e 


m 2 (say) 


(5.17) 


Here, all the field constants (B is ) are given at Appendix C. 



5.1.2 Numerical Implementation 


Using the same approach of numerical implementation as used in 
chapter 2, we can find radial wave number (k r ) and axial wave number (Pi) 
and thereafter resonance frequency can be computed. A software package 
(DRTM10) has been developed for this implementation and all the results 
generated in chapter 2, have been again verified, thereby concluding the cor- 
rectness of this analysis. Here also for top and bottom circular plates made 
of PEC, we have to use the modified DWM method as was done in section 
2.3. Similar to equation (2.33) we get 

( + i 7713 tan (M) _ tdjh ~ t&n(P 1 h) \ 

\e s .iPim 3 - e d p 2 tan(/Ah) e sl /3 1 m 4 + e d p 3 tan(/A/i ) ) ' 

( + €j4 ^ im3 tan (A^) _ tjsP 3 ~ CsiPim taxi{pih) \ _ 

\e si pim 3 - e ss p 2 tan (fah) e 5 i/3im 4 4- e ss p 3 tan^h) / - ' 

5.1.3 Other Derived Structures 

As it was mentioned earlier, using the structure of Fig 5.1, we can 
derive more structural environments for DR in TMois mode, out of which 
isolated DR, DR post, DR in MIC environment and DR in a cylindrical cavity 
have already been discussed in chapter 2. The other structures are given in 
Fig 5.2 to Fig 5.3, for which the resonance frequency can be computed using 
same analysis and software package. 


5.2 Energy Distribution 

The electrical and magnetic energy distributions in each region of Fig 
5.1, can be determined using the exactly same approach as done in chapter 3 
for six regions structure of Fig 2.1. Therefore, here we will present only the 
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Figure 5.3: 



5.2.1 Electrical Energy Distribution 


The electrical energies in each- region, are found to be as 
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Here constants B is are given at Appendix C and other constants are given 
at Appendix B. 



5.2.2 Magnetic Energy Distribution 


Similarly, the magnetic energy distribution in all the regions, is found 
to be as :- 
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5.2.3 Numerical Implementation 

Using the electrical and magnetic energies formulated in previous sec- 
tions, we can easily obtain the energy filling factors thereby finding the per- 
centage energy in each region, as done in chapter 3. Also, the total electrical 
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Figure 5.4: Edge coupled DRs placed in an MIC environment 


and magnetic energies can be computed and used in determination of Q- 
factor, inter-resonator coupling etc. 


5.3 Inter-Resonator Coupling 


Using the versatality of the structure in Fig 5.1, it is possible to de- 
termine the coupling between two adjascent DRs, kept in environments like 
in Fig 2.11(b&c), which can be called as broad-side coupling and in Fig 5.4 
which can be called as edge coupling. Skalicky [19] has reported the method 
for computation of inter-resonator coupling for TEW mode using EM in- 
duction and reaction concept. He has assumed the DRs operating in TE 0lS 
mode to be acting as magnetic dipoles (i.e. current carrying loops) and 
thereby induction of voltage into adjascent second loop, due to the current 
in the first one, using the inductive coupling between the two DRs. A de- 
tailed analysis based on this method has been carried out in [6] for TE 0lS 
mode, for determination of both broad-side and edge coupling between two 
adjascent DRs. Byzery [20] has used the similar approach for fidmg the edge 
coupling between the two adjascent DRs kept on a substrate, m rectangular 
waveguide (below cut-off) longitudanally, operating in TM 0 is mode. On the 
similar approach, a formulation has been done to find inter-resonator cou- 


plings between adjascent DRs kept in environments of Fig 2.11(b&c) and Fig 
5.4. The relation used in [20] based on [6] is given by :- 



(5.45) 


where k is the inter-resonator coupling, W e is the total electrical energy stored 
and E z \ & E z 2 are the longitudanal components of electric fields outside the 
first resonator and inside the second resonator. 


5.3.1 Edge Coupled DRs 


To find the inter-resonator edge-coupling between two adjascent DRs, 
shown in Fig 5.4 (a -4- 00 i.e. very large say « 5 d), we need to find the 
volume integral in equation (5.45). 



• E z2 dv 


i h n 

J-h JO Jo 


E z %E z i 


r + (2d -I- s) sin </> ' 


r dr d<f> dz 

(5.46) 


where, 


r' = |r 2 + (2d + s) 2 + 2 (2d + s)r sin <^j 1/2 (5-47) 

Here, s is the edge-to-edge inter-resonator gap, as shown in Fig 5.4. Here E z \ 
will correspond to the E z field in eight region of DR1 and E z2 will correspond 
to the E z field in the first region (DR) of DR2. On substituting the values 
of E sB and E zX from section 5.1.1 into equation (5.46), the volume integral 
can be evaluated as :- 

j ^ E zX ■ E z2 dv = B\qA\§Iko (5.48) 
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hi = /i 3 =0.4mra, lx -2 — /i 4 =l. 575mm, e 4 i = e s2 =2.54, e d = 24, Th=8.25mm 

Figure 5.5: Comparison of Inter-resonator edge-coupling with data published 
in [20] 

5.3.3 Numerical Implimentation 

A software package (DRTM10) has been developed for the entire 
analysis of the generalised structure of Fig 5.1. The results of inter-resonator 
coupling have been found to be very encouraging. Results of remaining 
parameters like resonance frequency, energy etc. are found to be the same 
as that in chapter2, and therefore those are not being shown here again. 


Edge Coupling 

The results for inter-resonator edge-coupling as a function of gap be- 
tween the resonators, shown in Fig 5.4, have been plotted in Fig 5.5 which 
have also been compared with the results published in [20]. 




is as expected because as substrate dielectric constant increases, more electri- 
cal eneigy gets stored in the substrate region and therefore more interaction 
with other DR takes place. Also, it is seen that Pure DWM method gives 
higher results for higher substrate dielectric constant, than that calculated by 
Present method, however the nature of variation of coupling with substrate 
gap, is same in both the cases. 

Fig 5.8 shows the variation in coupling between DRs through the 
air, as shown in Fig 2.11(b). This coupling has been calculated using both, 
Present method and the Pure DWM method. It is seen that this coupling 
however reduces with the increase in the substrate dielectric constant, be- 
cause as the substrate dielectric constant increases more electrical energy gets 
attracted towards the substrate region and lesser energy gets stored in the 
air-gap region and therefore less interaction with the other DR takes place. 
Here also, the Pure DWM mehod gives higher results than that obtained 
from Present method, however the nature of coupling variation is same in 
both the cases. 



Coupling (k) <S # Coupling (k) 



Substrate Gap Between DRs (mm) 


re 5.7: Variation of Inter-resonator coupling with substrate thickness 
/i 2 =1.0mm, /i 3 = /i 4 =5.0mm, e d =24, d=3.635mm, 
a=5.45mm, e s2 = e S 3 = e S 4=1.0, Th=4.04mm 



Air Gap Between DRs (mm) 

/?2= 1.0mm, /13 = /i4=5.0mm, e<j=24, d=3. 635mm, 
a=5.45mm, e s2 = e S 3 = e 5 i=1.0, Th=4.04mm 


Figure 5.8: Variation of Inter-resonator coupling with Air gap 




Chapter 6 


Filter Design 

C.l TM Mode DR Filter in Circular Wave- 
guide Below Cut-off 

Bandpass filters constructed by placing TEm mode DRs co-axially 
in a TEqi mode cut-off circular cylindrical waveguide, have been reported 
[21], For filters using the TM 0 is mode of DRs, on the other hand only few 
investigations have been presented [16], [20]. Filters using TM 0lS mode have 
been reported to have low loss operation [16] and a reduced excitation of 
spurious resonances as compared to TEqis mode filters [20]. In fact input 
and output couplings, realized by means of simple circular co-axial probes 
of suitable dimensions, allow excitation of only /'.A/o,, modes in the ciiculai 
waveguide, therefore avoiding spurious resonances of TEo n modes as well as 
those of hybrid modes for the circular symmetry of the structure. However, 
the choice of TM 0U resonance results in a longer structure than that obtained 
with the TE mi mode, for the same waveguide diameter. 



circular waveguide 



Input 

probe 


Output 

probe 


(b) 


Figure (5.1: Filter Structure (a) Cross-sectional view (b) Longitudanal view. 


6.2 Design Aspects 


The filter structure having axially symmetric DR placed in a circular 
waveguide below cut-off will have a cross sectional view as shown in Fig 6.1. 
The excited TM mode in the structure must be above cut-off inside the 
dielectric loaded sections and below cut-off in the empty waveguide sections, 
thus the constraint on waveguide radius a is, that it must be smaller than a 
maximum value a max given by 


vo 


hnax 


2 7T /() 


X 


(6.1) 


where, v Q is the velocity of light in the free space, /„ is the filter bandpass 
center frequency and x is the eigen value of the TM 01 in the empty waveguide 

(= 2.40482 ). 


at 



0.2.1 Criloriaii lor Selection of Parameters a and d 


1 lie relevant targets to be achieved by proper design are ■- 


1. Adequate limitation of the filter passband losses 

2. Adequate reduction of spurious passband 

3. Adequate reduction of overall filter structure size to achieve high degree 
of miniaturisation 


The criterion for selection of parameters a and d are based on compu- 
tation of unloaded quality factor (Q 0 ) and on spurious resonances. A good 
compromise between Q 0 , spurious resonances and structure overall size may 
be achieved, with presently available dielectric materials, by choosing a and 
d in the range of values (for = 20 - 40) [22]:- 


O.o Ojfiax ^ u ^ 0.6 Ojtjiolx 
1.1 1.3 

, — Umax ^ ® _ , — Umax 

V<7l y/td 


(6.2) 


6.3 Filter Synthesis 

Bandpass filter synthesis is usually started by designing a prototype 
lowpass lumped element network, using the insertion loss method. This 
lowpass prototype design has been the basis of design for a vast majority of 
filter [3, 23]. In insertion loss method, the design of the filter starts with 
the specifications about the insertion loss or the return loss for a lossless 
network, in the desired frequency band. After specifying the magnitude of 
the insertion loss as a function of frequency in the passband (Butterworth 
function, Chebvchcve function etc.) the prototype values are obtained from 
the formulas given in [23], Having obtained these prototype values the low 
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Figure C.2: Equivalent circuit representation of direct coupled bandpass filter. 

pass characteristic is transformed into bandpass characteristic by using well 
known lowpass to bandpass mapping. 

Bandpass filter which are physically symmetric can be realized as a 
simple combination of resonators. For such filter (excluding the generalised 
filter which uses cross coupling amongst resonators) Matthaei[23] has given 
an equivalent circuit shown in the Fig 6.2. Input and output couplings for this 
circuit is represented by external quality factors Q ex t(a) and Q e xt(b) while cou- 
pling between resonators is represented by coupling coefficient kjj+ 1 . These 
quantities are calculated in terms of low-pass prototype gu values from the 
formula given below 


Qexl(a) 

Qcxi(b) 

+ 1 


A r i 


(ffoV-R.) 

A'„ 


w 


9n9n - fl 

{Ki*+JR*) ~ w 




j = 1,2, 3...n- 1 




s/SjSj+i 


(6.3) 










here, w is t he fractional bandwidth given by 


W2 — U>i 
W — 

l(’o 

w»0 = V w 2 Wi (6 4) 

The passband edges w 2 and w u are for Chebycheve response, defined by pass- 
band iipple level. With Q ex t and k computed, design data are generated to 
relate these values to the electrical and physical parameters of the resonators 
and the inter-element spacing between them. This completes the synthesis 
of bandpass filter. 


6.4 Proposed Filter Design and it’s Specifi- 
cations 

Filter structure as shown in Fig 6.1 was concieved for design with 
the following specifications :- 

1. Response : Chebychev response with 0.5 dB ripple in passband 

2. Resonance frequency = 11.83 GHz 

3. 3 dB bandwidth = ±50 Mhz 

4. Stopband attenuation = 45 dB 

On choosing disc type resonator and selecting the parameters a and 
d as per the previous section, following electrical and physical parameters of 

Dlls are chosen 
Height 2 It — 4.04mm 
Diameter 2d = 7.27mm 
Dielectric constant of DR ed = 24.0 



Unloaded Quality factor Q 0 = 25000 

For these specifications the design data obtained is as follows :■ 
Number of Resonators = 4 


!h) 

= 1.0 

// 1 = 1.6703 

9 2 

= 1.1926 

(7s = 2.3661 

<7i 

= 0.8419 

g b = 1.9841 


Q ext (a) 

Q ext(b) 

9o9i 

U) 

9-i 95 

C J 

= 197.59 

= 197.6 


OJ 

= 5.989 x 10~ 3 

M/2 — 

VfjyJte 

A’2,3 = 

UJ 

= 5.032 x 10“ 3 

s/Filii 

A'3,4 = 

U) 

= 5.989 x 10" 3 



(6.5) 


The design charts for inter-resonator coupling for the DR specifica- 
tions chosen above, are given in Figs 6.3 for k^ 2 , k 3A and in Fig 6.4 for k 2 , 3 , 
keeping substrate thickness equal to 7.65mm. However the values of Q ex t a 
and Qextb will have to be found by measurement, as their formulation could 
not be done due to time constraint and complexity involved. 


n*7 
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Chapter 7 


Conclusion and Future 
Developments 


7.1 Conclusion 

The characteristics of Dielectric Resonator operating in TMq XS mode 
lias been studied in detail in various useful environments like Isolated cases, 
Post case, in cavity made of PEC/PMC circular plates, MIC enviroment, 
Suspended substrate environment and DR placed in a cylindrical waveguide 
below cutoff. The design parameters like resonance frequency, field con- 
figurations, energy distributions, computation of overall unloaded Q-factors, 
coupling between two resonators have been presented. The complete analysis 
is based on simple and approximate DWM method which gives reasonably 
accurate results, especially for DR post and Isolated DR. However, it was 
observed that pure DWM method gives high results for resonance frequency 
for close proximities of PEC surface, whereas for very close distances Mod- 
ified DWM approach is more accurate. As the distance of PEC surface is 
increased, EDC approach with little modification gives reasonable accurate 
results. The energy distribution and computation of unloaded Q-factor has 



boon done using basic formulations. The effect of various structural parame- 
ters on resonance frequency, energy distribution, Q-factor and inter-resonator 
coupling have been discussed in detail. Interestingly, the substrate dielectric 
constant which does not affect the resonance frequency so much, plays an 
important role in electric energy distribution in various regions of the struc- 
tiuc, and thereby it affects the quality factor and inter-resonator coupling 
tremendously. 

Ii/(h'( f ive user friendly software packages have been developed using 
the above analysis which can be used for CAD of microwave circuits like 
filters, oscillators etc. The generated results have been compared with the 
available published data. 


7.2 Future Developments 

The behaviour of resonance frequency for close proximities of PEC 
surfaces can further be analysed in order to obtain accurate formulation 
covering PEC plate distance from post case to isolated case. The input 
and output couplings are required to be formulated, and thereafter practical 
realisation of proposed filter structure can be carried out. TM 0 is mode in 
suspended substrate environment, has rarely been reported and therefore it’s 
applications in microwave circuit designs may be analysed further. 
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Appendix A 


Six Regions Structure 


A.l Values of Field Constants B is 


Aq 

MM 


A 

coth (cv a /i 2 ) 
tanh (a ft /?. 2 ) 

tan (A /i) -e d a a 

A 

coth (a a /i 2 ) 
tanh (a a /i 2 ) _ 

+ e<*a a tan (A A 


e d {cos (Pih) + Bi sin (A A} 


Q-Ckahi 

cosh (a a /i 2 ) 
sinh (a a /i 2 ) 


e d {cos (A/») - B\ sin (A ft)} 

e s e-“* /l3 

cosh (a, A) 
sinh (a,/ii) 


b 2 • b 6 

Jo (^fJ) 

/to (^aJ) ~ 

- B / 0 (Ad) 


B 7 = Z? 3 ■ B 6 


where, 


h.i = /) + I>2 


(A-l) 

(A.2) 

(A-3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 


(A.8) 
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/);, = ll + l>\ 

(A.9) 

A. 2 Ajjx and t Terms 


f h 

Aio = / {cos (ft z) 4- B\ sin (ftz)} 2 dz 

J « h 

(A. 10 ) 

B . * iilM (,- b8 ( 

(A.ll)’ 

r/i 

A U v = / {sin (ft z) - ftcos(ftz)} dz 

J - /i 

(A. 12) 

- ' i {' + - “%5 A) ( i - s o} 

(A.13) 

*(/.*,/ u,««) = e ” 2a “ /M / 4 shill 2 {a a (z - ft)} dz 

J h 

(A.14) 

sinh (2t* a /i 2 ) - 2a a /i 2 

4 a a e 2a °' u 

(A. 15) 

= e_2aa ' 14 cosh 2 {a a (2 — ft)} dz 

(A.16) 

sinh (2 a a /i 2 ) 4- 2a a /i 2 

(A. 17) 

4 a„ e 2a °' 14 

/(h lt h,,o.) = e~ 2o “ h3 [ sinh 2 {a s {z + h 3 )} dz 

J ~hz 

(A.18) 

sinh (2at s /ii) — 2a s ft 

4 a, e 2Q>/l3 

(A.19) 

r — h 

t' . = e~ 2a ’ h3 / cosh 2 {a^ (z + ft)} dz 

(A.20) 


sinh (2a s /ii) + 2 a s hi 
4 a 3 e 2a ‘ l13 
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Appendix B 


Integrals involving Bessel 

Functions 


(B-l) 


I JO ~ f r Jq ( k r r ) dr 

J o 

= - { ./? (Kd) + JS (Kd ) } 

//•;./ = [ r J'f (krr) dr 
Jo 

J { J'l (Kd) + 4 (Kd) - ~- d Jo (Kd) ■ Jr (Kd)) (B.2) 


! h in ~ J tl 1 { Ko(Kr) ~ BIo(K a r)} 


dr 


4 B 2 


■ {/o (k a r) - l\ ( k a r )} 


’-{K; (Kr) - I<l (Kt)} 

— B |r 2 {/o (k a r) Kir (k a r) + li (k a r) K\ (fc a r)}J 

Ikk, = I" r { A'i M + Bh (K.r) } 2 dr 

Jd 

’^{Kt(Kr) - Kl(K r) + ^K 0 (k a r) ■ Kr(Kr)) 
-{ 4(Kr ) - Ir(Kr) - ^Io(Kr) h(Kr)) 


d 

(B.3) 


i a 
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•f B 2 
4 B 


[ h M Ko (har) + j ^ 1 o M K ' M + h K] 


Ik 


o = /' / Jo (hr) {*, (kj) -BI„ for 1 )} sin ^ . 

JQ JO f 1 

(B.4) 


where, 


t — 4- (2d + S ) + 2 (2d + S) r sin <f> 


(B.5) 


This integral is evalvated using Nag subroutines (D01DAF). 


dr dcj) 



Appendix C 


Ten Regions Structure 


C.l Values of field Constants B is 


B 

Bi 

B 2 


B 3 


B \ 


B 5 
B h 
b 7 


Kq (fc a a) 

Iq ( k a a ) 

£d 02 + tan (Pih) 

e S 4 0i m 3 - e d 0 2 tan (0ih) 

01 { Bicos(Pih) — sin(/?i h) } 

02 { B 3 cos (0 2 h) - sin (0 2 h) } 
m x tan (0 2 h 5 ) — 1.0 


= -B 


mi + tan (0 2 h 5 ) 

02 { B 3 cos {0 2 h 5 ) - sin (0 2 h$)} 


a s3 e ~ a>3h6 


sinh {a s3 h 4 ) 
cosh (0:53/14) 

0i { Bi cos (0i h) + sin (0i h) } 

03 { Be cos (0 3 h) 4- sin (0 3 h) } 

1.0 — m 2 t&n(0 3 h 7 ) 
m 2 + tan (0 3 h 7 ) 

0 3 {Be cos {0 3 h 7 ) + sin {0 3 h 7 )} 


Be 


a s 2 e _Qj2/l8 


sinh (a s2 h 3 ) 
cosh (a s2 h 3 ) 


Bg = B 4 ■ Bio 
Bg — B 2 B\o 


(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 

(C.7) 

(C.8) 


78 


(C.9) 

(C.10) 



Jo ( k r d ) 

Ao ( k a d ) — B Iq ( k a d ) 
B u = B 5 , Bi 0 

B\2 = Bj ■ B\o 
where, 

/15 = h + /12 

/?6 = /l + /l2 + /I4 

h 7 =z h + hi 

h 8 = /i + /ii + /13 


C.2 and t Terms 


{cos (Piz) + Bism(Piz)} 2 dz 


sin (2Pih) 


= A { 1 + s; + ^ (l _^)} 
r/i 0 

= / {sin (A- 2 ) — Bi cos {Piz)} 2 dz 

- 4 + b ; - ^ (1 - *■)] 


/i 1 + B{ 


sin {2/3ih) 
2(hh 


(C.ll) 

(C-12) 

(C.13) 


(C.14) 

(C.15) 

(C-16) 

(C.17) 


(C-18) 

(C.19) 

(C.20) 

(C.21) 


An = [ 5 [B 2 {cos (faz) + B 3 sin (/? 2 z)}] ■ [B 2 {cos (&*) + B z sin (/3 2 jz)}]* dz 

J h 

(C.22) 

1 

(i) If B 2 , B 3 , p 2 are all Real, OR P 2 , B z Real and B 2 Imaginary then 

An - \Blh 2 [l + B\ + sm ^ 2/l2 l {(i _ B\) cos {P 2 (2 h + h 2 )} + 

2 P2 n 2 

2B Z sin {P 2 (2/i 4- /i 2 )}] (0-23) 
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(ii) If 13 2, B 3 Real, and @2 Imaginary OR B 3 Real and /? 2 , B2 Imaginary 
then 

= \ B lH 1 ~ B l + 2) {C 1 + B D cosh {fe (2 h + hi)}}] 

(C.24) 

(iii) If p 2 , B2 Real, and B$ Imaginary OR p 2 Real and B2, B 3 Imaginary 
then 

-1,7 = \B 2 2 h 2 [l + Bl + S -^^-{(l-Bl)cos{l3 2 (2h + h 2 )}}tG.25) 

(iv) If B 2 Real, and fe, B z Imaginary OR /? 2 , B 2 , B 3 are all Imaginary 
then 


.4 17 — -Blh 2 [l — Bl 4- — -^-{(1 + B\) cosh {P 2 (2/i + h 2 )} — 

2R3 sinh {/?2 (2/r + /12)}] (C.26) 


A 'n = [ \B%fe [B 3 cos (fez) - sin (fez)}] ■ [B 2 p 2 [B 3 cos (fez) - sin (fez)}]* dz 

J h 

(C.27) 

(i) If B 2 , B 3 , p 2 are all Real, OR p 2 , B 3 Real and B 2 Imaginary then 

A'„ = \BlPlh 2 [l + Bl - 1 - Bl) cos {A (2 h + A)} + 

2B3 sin {/?2 (2/i + /12)}] (C.28) 

(ii) If B 2 , B 3 Real, and fe Imaginary OR B 3 Real and fe, B 2 Imaginary 
then 


,4 


/ 

17 


^BlPlh 2 [Bl - 1 + SinI ^" ~{( 1 + B D cosh ( 2h + 

(C.29) 
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(iii) If 02, B 2 Real, and B 3 Imaginary OR 0 2 Real and B 2 , B 3 Imaginary 
then 

j4 i7 = + B\ + — ^^3 ~ !) cos {02 (2 h + /z 2 )}}] 

(C.30) 

(iv) If B 2 Real, and 0 2 , B 3 Imaginary OR 0 2 , B 2 , B 3 are all Imaginary 
then 

A\ 7 = -Bl0lli 2 [Bl - 1 + — — {(1 + B]) cosh {0 2 (2 h + h 2 )} - 

2 B 3 sinh {0 2 (2 h + h 2 ) }] (C.31) 

Note Here only neumerical values of real or imaginary parts of B 2 , B 3 and 
0 2 are to be substituted. 

r~ h 

-^18 = / [Bs {cos (0 3 z) + B 6 sin (0 3 z)}]-[B 5 {cos (0 3 z) + B 6 sin (0 3 z)}]* dz 

^i8 = f [B 5 P 3 {R6 cos ( 0 3 z ) - sin (0 3 z)}} ■ [B 5 0 3 {B 6 cos ( 0 3 z ) - sin (0 3 z)}]* dz 

J —hi 

(C.32) 

Values of these integrals can be obtained by replacing B 2 by B 5 , B 3 by jB 6 , 

P 2 by 03 and h 2 by h 3 in the values of A u and A' u with signs of last two 
terms reversed in case (i) and (iv). 


Hh<,h 6 , a , 3 ) = e 2033,16 / sinh 2 {a S 3 (2 — /i 6 )} dz (C. 33 ) 
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4 a s3 e 2a ’ zh 6 
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J — h& 


81 



f r 

b (hz,ft&,Oi s 2) 


(C.38) 

(C.39) 

(C.40) 


sinh (2a s2 /i 3 ) - 2tt, 2 ft a 
4 a 52 e 2a * 2h & 

e ~ 2a,2h * f_ ha cosh 2 {a s2 (z + h s )} dz 
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4 q s2 e 2a ° 2h 8 
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= - B X sin( M } • sinh{a 42 (ft 3 - ft + z)} dz 

= /?! + a£> ^ ^ Sin ^^ {si n h(a s2 / l3 )+sinhK 2 (h3-2/z)} } + 
a s2 cos(/?i ft) {cosh (a s2 ft 3 ) - cosh {a i2 (ft 3 - 2h)}} - 
B 1 {A cos (Aft) {sinh{a s2 (ft 3 — 2h)} - sinh(a s2 ft 3 } + 
a s2 sin{Pih) { cosh {a s2 (ft 3 - 2ft)} + cosh(a j2 ft 3 } }] (C.41) 
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